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Abstract — The compound models of clutter statistics are found 
suitable to describe the nonstationary nature of radar bacltscat- 
tering from high-resolution observations. In this letter, we show 
that the properties of Mellin transform can be utilized to generate 
higher order moments of simple and compound models of clutter 
statistics in a compact manner. 

Index Terms — Clutter, compound model, Mellin transform, 
log-cumulants. 



I. Introduction 

RADAR backscattering from ground or sea surfaces are 
wide-sense stationary for low-resolution observations as 
expectations of clutter statistics or moments are assumed to be 
independent of spatio-temporal changes. For high-resolution 
observations, such surfaces reveal heterogeneous structures 
such as swell in sea waves or winds blowing over the canopy 
of grasslands that result in nonstationary clutter statistics [1], 
[2], [4]. The compound models of probability density functions 
(pdf) incorporate the variation in the parameters of clutter in 
such cases. 

Traditionally higher order moments of a continuous random 
variable (rv) X are generated from higher order derivatives of 
its characteristic function defined as 



^xit^) = E{e^p(jujx)} 



fx {x) exp (jujx) dx. ( 1 ) 



The continuous pdf^(x) is for — oo < x < oo. Generation of 
moments and cumulants from ([T]) for the compound models of 
clutter require solutions of incomplete integrals. The domain 
of Z is < jc < oo for amplitude and power statistics, 
and /p fx{x)dx — 1. Properties of Mellin transform provide 
the formalism to derive higher order moments in a compact 
manner in such cases. Some of these properties were used in 
[3], [6] to derive the moments for high-resolution synthetic 
aperture radar (SAR) clutter statistics. Here, we show that the 
properties of Mellin transform can be utilized in an effective 
manner for both simple and compound models of clutter either 
in amplitude or in intensity domain. 

II. Mellin Transform Properties 

Mellin transform exists for a continuous function fx{x) 
defined over R_|_. The transform operator is the second kind 
characteristic function $^"(5) expressed as 



/■oo 

^xis)=M\fxix);s]= / x^~%{x)dx. 
Jo 



(2) 



Here s = a+jb G C is the complex Laplace transform variable. 
Traditional moments are generated from ^ with s = n+l,n E 



m„ = M 



(3) 



s=n+l 



Second-kind moments or the log-moments are generated for 
logarithm of rv X by using the derivative property of Mellin 
transform. 



in„ = M[log(x)%{x);s] 



x'-^\og{x)%(x)dx 

s=l "'0 

^M^) . (4) 



Analogous to the cumulants derived from logarithm of 
characteristic function in ([T]), the n-th order cumulants of 
second kind or the log-cumulants are obtained from derivatives 
of logarithm of ^x{s); i.e., *a-(s) = log((f>A-(s)). 



(5) 



The log-moments and the log-cumulants are related as 



ki = rill 



k2 — ni2 — mi^ 
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— m^ — 3mitn2 + 2mi 
k^ — TO4 — Arfiim^ + 6mi^n'i2 — Snii'^. 



(6) 



The underlying mean of speckle component of clutter 
vary widely in the compound models of amplitude or power 
statistics resulting in long-tailed distributions. Speckle arises 
from randomness in the distribution of backscattering elements 
in the resolution cell, the number of such scatterers is non- 
stationary for high-resolution observations. The pdf of high- 
resolution clutter is described by taking into account of a rv 
Z signifying randomness in the mean of clutter 



fxix) 



fxix\z)fz{z)dz, 



z>0. 



(7) 
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The compound pdf model in (Q is a Mellin convolution. One 
nice property of Mellin transform is the product form of the 
components of pdf in the transform domain [5]. 



(8) 



The log-cumulants of the components in ([Sj are therefore 
additive. 



(9) 



III. Moments Generation for Simple Models of 
Clutter 

The shape and scale parameters of simple models of pdf for 
low-resolution cases are stationary. The usual pdf of speckle 
power is a gamma distribution resulting from convolution of 
L independent exponential distributions. 



/v(v) 



1 



r(L)(/x 



v,(i-i)exp( ), 



V > 0. (10) 



Here r(.) is the standard gamma function. The shape and scale 
of distribution are determined by L and /i, mean value of clut- 
ter power respectively. Corresponding amplitude distribution 
turns out to be a Nakagami pdf [2], [6]. 



exp 



r(L) V 

Mellin transform for gamma pdf is 



LP 



r>0. (11) 



r(i) Jo 



(L+j-l)-l 



exp(— Av)Jv (12) 



with A = -. Using the transform pair 



M[x" exp(-Aji:);i] 



!-(■'+") 



r(i + u) 



we obtain. 



We use the additional Mellin transform pair 



1 



X[cxp(-Ax^);.]^-(A)-^r(-) 



with A — 2^ ; so that 



The moments of first kind for Maxwell pdf are 



(16) 



1 /2 



-W-(2a^)^ 

(7 V TT 



n + l\ ^/n + l 



(17) 

The moments for amplitude distributions are also derived by 
Mellin transform. As for Nakagami distribution, with A = — 



^^-(.-i)r(L + £_i) 



exp(— A^r^)t/r 



r(L) 



(18) 



The log-cumulants are easier to derive here. In general the 
log-cumulants of Nakagami distribution are derived from (|5]l 



T(n- 1,L) 



(19) 



Here T(.) is the Digamma function; i. e., the first derivative 
of In r(i) at i = 1. In general T(n — 1, L) is the nth derivative 
of the Digamma function for variable L. 

One long-tailed pdf often used in sea-clutter amplitude 
modelling [1] is Weibull distribution. 



b\ X\ r X 



cxp [-( -J > ^;b,z > 0. 

(20) 



Here z is the scale parameter and b is the shape parameter of 
distribution. Mellin transform of (l20ll is 



(13) 



The moments of first kind for gamma pdf are generated from 
(fTsT i with s = n + 1 as 



j=«+i 



m ■ 



(14) 



As a special case of the result in (fl4] l. the moments of 
exponential pdf (for L — 1) are m„ = fj."n\. Maxwell pdf 
is the case for L = 3. 



/M(«) = ^y§«'exp(-^), «>0. (15) 



$^(,) = ^^°°x('+''-i)-iexp[-Q ]dx. 
From the Mellin transform pair 

b 

the second characteristic function is 

,s + b-l. 



The moments of first kind for Weibull distribution are 

„ /n + b^ 

, , b 



(21) 



(22) 



(23) 
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The common Rayleigh amplitude pdf is a special case of 
WeibuU distribution with 6 — 2. 



/«(r;z) = 2(^^)exp[-Q ], 



r > 0. (24) 



The moments of first kind for Rayleigh pdf are 

'n + 2 



m„ = $if(s) 



z"T 



s—n-\-l 



(25) 



is a Mellin convolution of Rayleigh pdf and exponential pdf 
given by. 



fNir) 



Art' 



Q pOG 



r(a) Jo 2 



z exp — Dz 



(30) 



Here a is the shape parameter in the variation of mean of 
sea or ground clutter amplitude, and b is the scale parameter 
for associated speckle amplitude. Following derivation for 
Nakagami pdf in (fT9] l the second characteristics function for 
K-pdf is given by. 



We show in the next section the utility of Mellin transform 
for deriving the log-moments and the log-cumulants of com- 
pound models of clutter in a compact manner. 

IV. Moments Generation for Compound Models of 
Clutter 

The pdf for compound models of high-resolution clutter 
have got two components; pdf of speckle component, and pdf 
of the modulation in mean amplitude or power of speckle. 
Considering both to be gamma distributed rv the pdf for 
generalized gamma (GF) model of clutter power is [6] 



./v(v) 



r(L)r(M) V< z > 



2LM 



2LM 



<z> 



Km-l [2 



LM 

<z> 



(26) 



The shape parameter for gamma pdf /2(z) of rv Z according 
to (I?) is M, and Km-l{-) is the second kind modified Bessel 
function of order (M — L). The mean estimate of < z >= /i. 
Assuming speckle and the modulation in mean power in the 
high-resolution cell to be independent of each other, we have 
by Mellin convolution property in ^ 



^y{s) = [s - 1) log j + logF(. + L - 1) + 

logF(i+M- 1) -logr(L)-logF(M). (27) 



The log-cumulants of GF model are 



^~i=log(^)+T(L) + T(M) 
k„ = T(n - 1,L) + T(n - l,M). 



s=l 



(28) 



Spikes in high-resolution ground-clutter amplitude at low 
grazing angles are often described by the K-distribution model 
[2], [4]. The compound K-pdf 



(c + l) 

4/7 ^~r" r- 

T{a) 



(29) 



where < z >— /x^. The log-cumulants for K-pdf are 



(31) 



and 



^^log^+(i-l)logAi + logF(^^ + i^ + 

logF(a), 



logF a 



s- 1 



ki = -^logb + logn+^T{a) 
4= T(«-l,a). 



-T(l) 



(32) 



Here T(l) = —0.577215 is the Euler constant [5]. This shows 
that the log-cumulants of K-distribution are determined by the 
higher order log-cumulants of Nakagami distribution in the 
mean of high-resolution ground or sea clutter 

A more extended case of compound clutter model is the 
scene where variation in the shape of clutter amplitude distri- 
bution is given by generalized WeibuU distribution [4], 



fwNinb,c,a) 



2cfe" 
Via)' 



.e-l 



2a-l 

2 



-bz^ 



(33) 

This is a Mellin convolution where randomness in the mean 
amplitude of clutter is described by Nakagami pdf with the 
shape parameter being a. 



2fe" 



,2a-l 



exp(-fe2), 



and the clutter amplitude follows a generalized Weibull distri- 
bution with the shape parameter being c. 



/H'(''k;c) = -r' 'exp 



Following the transform rule of Mellin convolution. 



F 



- 1\ Fa 



F(a) 



(34) 
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where < >= a. The log-moments for this generalized 
WeibuU model of clutter according to (|9]l are 



1 / a 
^1=2 



-Y(l) + -Y(a) 



K=\\ \ T(n-l,a). 



(35) 



Another compound model used to describe high-resolution 
SAR clutter is the Fisher distribution [3]. 



/f(«) = 



r(L + M) / L \ (sItm) 



r(L)r(M) Vm/.; + «) 



, L+M- 



(36) 



Consider = A; following the Mellin transform pair 



x[(l+A)-^. 



with b ^ L + M 



the second characteristic function for Fisher distribution is 



E ' 

o 

E 
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Fig. 1. Log- moments of data and log-cumulants of texture for simulation 
of GF model of high-resolution clutter. Left: second order moments; right: 
fourth order moments. 



1 



V{s + L-l)V{M+l~s). (37) 



r(L)(M) 
Corresponding log-cumulants are 

ki log^i + [T(L) - logL] + [Y(M) - logM] 
k„ = T(n - 1,L) + (-l)"T(n - 1,M). 



(38) 



One useful application of the log-cumulants and their 
relationship with the log-moments in (|6]l is estimation of 
parameters of texture. Empirical data from high-resolution 
radar backscattering x{t) follow the product model, 



x(t) = u{t)z{t). 



(39) 



Here u{t) is the speckle component, and z{t) represent texture 
signifying variation in the mean parameter The log-moments 
of observed data and the log-cumulants of texture can be 
estimated for different compound models utilizing the relation- 
ships in (IHi and (|6]l. Parameters for texture are derived using 
the log-cumulants of speckle as in ©, and can be verified with 
the theoretical values derived in the paper For example, second 
and fourth order log-cumulants of texture component for GF 
model of high-resolution ground clutter in ( l28T l are estimated 
in Fig. 1. Second and fourth order log-moments of x{t) are 
derived from the log-cumulants of z{t) assuming it to be a 
gamma variable, and u{t) also follows gamma distribution. The 
results of simulation show that higher order log-cumulants of 
texture vanish with increasing values of shape parameter M. 
This is expected in the present case as the texture component 
follows a nearly Gaussian distribution with constant mean 
for increasing values of M. For values of M < 1, there is 
presence of large amount of spikes in observed data signifying 
high values of log-moments and cumulants. Log-moments of 



clutter tend to become constant with increasing M signifying 
stationarity of low-resolution observations. 

V. Conclusion 

The utility of Mellin transform properties to generate 
higher order moments of simple and compound models of 
clutter in both amplitude and power domain is shown in 
this letter. The second kind characteristic function and its 
properties provide compact analytical expressions for higher 
order moments that are useful to interpret texture properties 
of high-resolution clutter. 
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